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ABSTRACT 

We present Goddard High-Resolution Spectrograph observations of the interstellar H i and D i Lya lines 
and the Mg ii and Fe ii resonance lines formed along the lines of sight toward the nearby stars Procyon (3.5 
pc, / = 214°, b = 13°) and Capella (12.5 pc, / = 163°, b = 5°). New observations of Capella were obtained at 
orbital phase 0.80, when the radial velocities of the intrinsic Lya emission lines of each star were nearly 
reversed from those of the previous observations at phase 0,26 (analyzed by Linsky et al.). Since the intrinsic 
Lya line of the Capella system— the “continuum” against which the interstellar absorption is measured— has 
different shapes at phases 0.26 and 0.80, we can derive both the intrinsic stellar profiles and the interstellar 
absorption lines more precisely by jointly analyzing the two data sets. We derive interstellar parameters from 
the simultaneous analysis of the two data sets as follows: (D/H)hsm = (^-^O ± 0.09 [+0.05, -0.10 systematic 
error]) x 10" ^ temperature T = 7000 ± 500 [±400 systematic error] K, and microturbulence = 1.6 ± 0.4 
[and ±0.2 systematic error] km s"^ (All random errors determined in this paper are ±2 (T.) 

For the analysis of the Procyon line of sight, we first assumed that the intrinsic Lya line profile is a 
broadened solar profile, but this assumption does not lead to a good fit to the observed D i line profile for 
any value of D/H. We then assumed that (D/H)l,sm = 1-6 x 10 the same value as for the Capella line of 
sight, and we modified the broadened solar profile to achieve agreement between the simulated and observed 
line profiles. The resulting asymmetric intrinsic stellar line profile is consistent with the shapes of the scaled 
Mg II line profiles. We believe therefore that the Procyon data are consistent with (D/H)lism = L6 x 10 but 
the uncertainty in the intrinsic Lya emission-line profile does not permit us to conclude that the D/H ratio is 
constant in the local interstellar medium (LISM). The temperature and turbulence in the Procyon line of sight 
are T = 6900 ± 80 (±300 systematic error) K and = 1.21 ± 0.27 km s“^ These properties are similar to 
those of Capella, except that the gas toward Procyon is divided into two velocity components separated by 
2.6 km s"‘ and the Procyon line of sight has a mean neutral hydrogen density that is a factor of 2.4 larger 
than that of the Capella line of sight. This suggests that the first 5.3 pc along the Capella line of sight lies 
within the local cloud and the remaining 7.2 pc lies in the hot gas surrounding the local cloud. 

We propose that = 0.1065 ± 0.0028 cm“^ be adopted for the neutral hydrogen density within the local 
cloud and that = 1,21 ±0.27 km s“‘ be adopted for the nonthermal motions. The existence of different 
second velocity components toward the nearby stars Procyon and Sirius provides the first glimpse of a turbu- 
lent cloudlet boundary layer between the local cloud and the surrounding hot interstellar gas. We speculate 
that what is often called “ turbulence ” may instead be velocity shear within the local cloud that is not a rigid 
comoving structure. We also derive gas phase abundances of iron and magnesium in the Procyon line of sight 
and the abundance of oxygen in the Capella line of sight. 

Within the context of standard big bang nucleosynthesis, our observed value of (D/H)l,sm leads to 0.042 ^ 
Db ^50 ^ 0.09, depending on the assumed model for Galactic chemical evolution of deuterium. Our lower limit 
(D/H)lism > L41 X 10“® provides a hard lower limit to the primordial D abundance and thus a hard upper 
limit on CIbHo ^ 0.125. These limits are independent of Galactic chemical evolution models and only assume 
that D is destroyed with time. 

Subject headings: ISM: abundances — stars: individual (a Canis Minoris, a Aurigae) — ultraviolet: stars 


1. INTRODUCTION 

The Hubble expansion, microwave background, and light- 
element abundances are the main observational pillars upon 
which the standard big bang cosmology now rests. Of the 
three, the light-element abundances provide the main con- 
straint on the total baryon density (both luminous and dark 
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baryonic matter). The D/H ratio is the most stringent of these 
constraints because (1) there are no known significant sources 
of deuterium after about 10^ s in the early universe, (2) the 
subsequent destruction of D is only by nuclear reactions in the 
cores of stars where D is the most fragile species, and (3) there 
is a steep monotonic slope between the primordial D/H ratio 
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and the baryonic density in contemporary big bang nucleo- 
synthesis models (e.g.. Walker et al. 1991). Since none of the 
other light elements (^He, "^He, ^Li, ^Li, Be, or B) share these 
properties, their abundances provide less certain estimates of 
the baryon density of the universe. 

The importance of measuring D/H has led to searches in a 
wide variety of environments including terrestrial seawater 
(using HDO/HjO), atmospheres of the giant planets (using 
HD/H 2 and NH2D/NH3), cold interstellar clouds (using deu- 
terated molecules and the as-yet-unsuccessful search for the 92 
cm hyperfine transition of atomic D), and warm interstellar gas 
(using the Lyman series lines of D i and H i). This work has 
been reviewed by Boesgaard & Steigman (1985) and most 
recently by Wilson & Rood (1994). Deuterated molecules have 
been observed in many sources, but their abundances relative 
to the corresponding undeuterated molecules are generally 
larger (often by orders of magnitude) than expected for sensible 
values of D/H. The overabundance usually is explained by the 
large mass difference between D and H, the slightly larger 
binding energy of the deuterated molcules, and the cold 
environments in which they formed (e.g., Geiss & Reeves 1981). 

Determining the ratio of the column density of D i to that of 
H I in warm interstellar gas (T % 7000 K), by measuring the 
absorption in the Lyman series lines along the lines of sight 
toward stars, is now thought to be the most accurate method 
for inferring the D/H ratio in the Galaxy. Although the inter- 
stellar gas has been chemically processed over the lifetime of 
the Galaxy and the D/H ratio must therefore be lower than 
primordial, this method should not suffer from other system- 
atic errors. In particular, the relative ionization fraction, 
molecular association fraction, and degree of condensation 
onto dust grains should be the same for D and H in this 
environment. For Galactic lines of sight, only the Lya line can 
be studied by the International Ultraviolet Explorer (lUE) and 
the Hubble Space Telescope (HST), but the higher. Lyman 
lines were observed toward a few hot stars by Copernicus. 
Unfortunately, H i absorption overlaps the D i line (—0.3307 
A or —81.55 km s“ ^ from the H i Lya line). This overlap limits 
the use of the Lya line to relatively nearby stars where the H i 
lines is not too opaque at the center of the D i line (log Nh i < 
18.7). A reanalysis of the best available Copernicus and lUE 
data led McCullough (1992) to estimate that the mean value of 
D/H by number in the local interstellar medium (LISM) is 


(D/H)l,sm = L5 ± 0.2 X 10"^. Unfortunately, these data have 
a relatively low signal-to-noise ratio (S/N), the H line is very 
saturated, and the intrinsic shapes of the D i and H i lines are 
unresolved at the spectral resolutions of lUE and Copernicus. 
We therefore initiated an HST observing program to obtain 
accurate values of (D/H)lism, to estimate the destruction of D 
over the lifetime of the Galaxy and to infer the primordial D/H 
ratio. 

Since 1990 the Goddard High Resolution Spectrograph 
(GHRS) on the Hubble Space Telescope has been acquiring 
ultraviolet spectra with unprecedented spectral resolution and 
S/N. (See Brandt et al. 1994; Heap et al. 1995; and Soderblom 
1993 for a description of the GHRS and its in-orbit capabilities 
at the time of these observations, prior to the installation of 
COSTAR.) The high spectral resolution and S/N of these data 
permit us to measure very accurate column densities of inter- 
stellar D I, H I, and the heavier elements, together with the 
thermal and dynamic properties of interstellar gas, in the 
region of space near the Sun that includes the local warm cloud 
and surrounding hot substrate plasma and additional warm 
clouds. (See Vidal-Madjar et al. 1978; Lallement & Berlin 
1992; Frisch 1994a, b; Linsky et al. 1993 [hereafter Paper I]; 
and below for a discussion of the LISM properties.) In this 
second paper in a series, we reanalyze the Capella line of sight 
in § 2 using observations obtained at both orbital quadratures, 
and in § 3 we analyze observations of the Procyon line of sight. 
We discuss our results in § 4 and present our conclusions 
in§ 5. 

2. THE CAPELLA LINE OF SIGHT 
2, 1 . Previous Observations 

Our program to measure the D/H ratio and the gas proper- 
ties in the LISM began on 1991 April 15 with echelle observ- 
ations of the resonance lines of H i and D i (Lya at 1216 A), the 
Fe II multiplet UVl (at 2599 A), and the Mg ii and k lines (at 
2796 and 2803 A) toward the Capella binary system. These and 
the subsequent GHRS observations of Capella (and Procyon) 
are summarized in Table 1. All the observations were obtained 
through the small science aperture (SSA) of the GHRS using 
the substep patterns listed in Table 1 (see Soderblom 1993; 
Paper I). The first observations were made at orbital phase 
0.26, when the G8 III star in the Capella system had a radial 


TABLE 1 

Summary of GHRS Observations 


Aperture Spectral Spectral Exposure Start 

and Substep Range Resolution Time Time Important 


Grating Pattern (A) (km s b (s) (CT) Features 

Capella at Phase 0.26 (1991 Apr 15) 

EA-46 SSA 9 1211-1217 3.57 3917 10:04 H M- D i >U2I6 

EB-22 SSA 7 2594-2606 3.27 707 13:37 Fe ii >12599 

EB-20 SSA 7 2793-2807 3.54 707 13:23 Mg ii A;12796, 2803 

Capella at Phase 0.80 (1993 Sep 19) 

G160M SSA 5 1206-1243 15 979 11:11 H i + D i >U216 

EB-20 SSA 7 2792-2807 3.54 707 12:48 Mg ii >1^2796, 2803 

Procyon (1992 Dec 21) 

GI60M SSA 5 1206-1243 15 3264 17:31 H i -I- D i >11216 

EB-20 SSA 7 2793-2808 3.54 653 19:22 Mg u XX2796 . 2803 

EB-22 SSA 7 2594-2606 3.27 653 19:47 Fe ii >12599 
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1215.0 1215.5 1216.0 1216.5 

Wavelength (A) 

Fig. 1. — Lya profiles of Capelia at 0 = 0.26 and 0.80. Irregular solid line, high-dispersion echelle-A spectrum obtained 1991 April 15 at phase (f> = 0.26. The 
profile consists of a broad Lya emission line produced in the chromospheres of both stars and absorption by interstellar H i and D i (A/i = -0.33 A relative to 
hydrogen) along the 12.5 pc line of sight to Capelia. Dashed line, 0 = 0.26 spectrum smoothed to the moderate dispersion G160M resolution. Smooth solid line, 
moderate-dispersion G160M spectrum obtained 1993 September 19 at 0 = 0.80. The difference between the smoothed 0 = 0.26 spectrum and the 0 = 0.80 spectrum 
outside of the saturated cores of the H i and D i lines is caused by the different radial velocities of the stellar emission lines at these two phases. 


velocity of +55.6 km and the G1 III star had a radial 
velocity of +2.0 km very close to their maximum 

radial velocity separation.^ The broad H i interstellar absorp- 
tion and narrower D i absorption are superposed upon 
Capella’s chromospheric H i Lya emission line (see Fig. 1). 
Detailed analysis of this spectrum and the Fe ii and Mg ii 
spectra in Paper I permitted us to obtain very accurate mea- 
surements of both the temperature {T = 7000 + 200 K) and 
nonthermal broadening = 1.66 + 0.03 km s"^) that charac- 
terize the gas in this line of sight and may be representative of 
the gas parameters in the local cloud. A careful analysis of the 
Capelia spectra, including systematic errors associated with 
the uncertain intrinsic emission-line profile against which the 
interstellar absorption is measured, shows that the neutral 
hydrogen column density is = 1.7-2, 1 x 10‘® cm“^ and 
the number density ratio (D/H)hsm = 1.65 ( + 0.07, 

— 0.18) X 10"^ for this line of sight. 

2.2. New Observations of Capelia at Phase 0.80 

The largest identified systematic error in our previous 
analysis of the Capelia phase 0.26 observations was the uncer- 
tain intrinsic Lya emission-line profiles of the two stars in the 
Capelia system, especially those parts of the emission lines 
within 0.5 A of line center but outside the dark core of the 
interstellar H i Lya absorption line. These parts of the intrinsic 

^ We use the ephemeris of Ayres (1984, 1988), in which phase 0.0 is in 
conjunction with the slightly more massive G8 HI star in front. 


line profiles are critically important because they form the 
“continuum” against which the observed absorption-line 
profile is compared to determine the interstellar column den- 
sities and broadening parameters for H i and D i. It was diffi- 
cult to solve for both the shape of the two intrinsic stellar Lya 
emission lines and the interstellar gas properties with spectra 
obtained at only one orbital phase, although lUE spectra at 
other phases were helpful (Ayres et al. 1993). We therefore 
reobserved Capelia on 1993 September 19 at orbital phase 
0.80, close to the opposite orbital quadrature (phase 0.75). Our 
goal was to verify the results in Paper I by analyzing the 
(assumed constant) interstellar absorption against the back- 
ground of a somewhat different intrinsic emission line from the 
Capelia system. At phase 0.80 the radial velocities of the G8 III 
and G1 III stars were +4.7 and +55.5 km s"*, respectively, 
nearly the reverse of their velocities at phase 0.26. Our observ- 
ations (see Table 1) consisted of a G160M spectrum of the Lya 
region (see Fig. 1) and an echelle spectrum of the Mg ii region 
(see Fig. 8 in Linsky et al. 1995). The echelle-A grating on side 1 
of the GHRS was not operational at the time of these obser- 
vations. These observations were obtained in the same way as 
described in Paper I, except that only the G160M spectrum 
was obtained with the FP-SPLIT procedure to reduce the 
photocathode fixed pattern noise and granularity. The on- 
board Doppler compensation to correct for the changing 
spacecraft radial velocity was enabled during all observations. 

The processing of the phase 0.80 spectra and reprocessing of 
the phase 0.26 spectra employed the 1993 December version of 
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the CALHRS calibration software at the GHRS computing 
facility at the Goddard Space Flight Center. The processing 
includes converting raw counts to count rates and correcting 
for paired pulse events and diode-to-diode nonuniformities. 
The quarter diode stepp>ed spectra were merged to form the 
individual spectra of 2000 wavelength samples. The final 
reduction includes proper alignment of the individual FP- 
SPLIT spectra (using the HRS MERGE procedure), rejection 
of fixed pattern noise structures, correction for scattered light 
and the radiation background, and wavelength calibration 
including the effects of thermal drifts. For the echelle-B mea- 
surements, the background correction was determined as 
described in Paper I. Wavelength calibration using the 
CALHRS reduction software as of 1993 December was per- 
formed with the special platinum lamp calibration exposures 
obtained immediately before each stellar integration. 

In Figure 1 we compare the observed Lya line at phase 0.80 
with the phase 0.26 profile degraded to the resolution of the 
G160M grating. The two profiles are nearly identical in the 
cores of the H i and D i lines where the interstellar optical 
depth is large and the shape of the underlying stellar emission 
line is relatively unimportant, but the two profiles differ greatly 
where the interstellar opacity is small, including spectral 
regions on both sides of the core of the D i line. The avail- 
ability of data at both quadratures makes it feasible to deter- 
mine both the shapes of the stellar emission-line profiles and 
the interstellar absorption parameters. 


2.3. Combined Analysis of the Phase 026 and 0.80 
Capella Spectra 

Our analysis of both Capella data sets differs from that used 
in the previous analysis of the phase 0,26 data in Paper I in two 
important respects. First, we developed models for the intrinsic 
Lya emission-line profiles of the G1 III and G8 III stars (see 
Fig. 2) and required that they be invariant with phase, aside 
from the different orbital velocity shifts. This requirement is 
sensible, since Ayres (1984) and Wood & Ayres (1995) have 
detected no significant changes in the Capella ultraviolet 
emission-line fluxes with orbital phase or time at the ±5% 
sensitivity level of IV E, During the iterative fitting process, we 
modified the interstellar parameters and the intrinsic line pro- 
files for each star, although the changes from the intrinsic pro- 
files used in Paper I are small, to minimize the value of 
between the observed and simulated Lya line profiles at both 
phases. 

A second difference from our original analysis in Paper I is 
the use of the two-component instrumental point spread func- 
tion (PSF) for the echelle data described by Spitzer & Fitz- 
patrick (1993). It consists of a narrow Gaussian “core” with 
FWHM = 1.05 diodes and a broader Gaussian “halo” with 
FWHM = 5.0 diodes. The relative strengths of these two com- 
ponents depend on wavelength (Cardelli, Ebbets, & Savage 
1990). We do not explicitly include the third component of the 
PSF measured by Cardelli et al. — the broad power-law 



1214.0 1214.5 1215.0 1215.5 1216.0 1216.5 1217.0 1217.5 

Wavelength (A) 

Fig. 2. — Capella at (^ = 0.26 and 0.80. Upper graph : High-dispersion (A/AA = 84,000) spectrum (displaced upward oy two units) of the 1216 A re^on of Capella at 
<p - 0.26 obtained with the echelle-A grating. Dashed lines are the assumed intrinsic emission-line profiles of the two stars in the Capella system. Thin solid line is the 
sum of the two intrinsic profiles folded through the interstellar medium with the parameters listed in Table 2. Lower graph: Moderate-dispersion {X/AX = 20.000) 
spectrum of the 1216 A region of Capella at </> = 0.80 obtained with the G160M grating. Dashed lines are the same intrinsic line profiles, but with the radial velocities 
of each star at this phase. 
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“wing” — but instead, following Spitzer & Fitzpatrick, we sub- We next reanalyzed the phase 0.26 Lya profiles by using the 

tract the measured flux in the saturated core of the interstellar two-component PSF described above. The observed and simu- 

Lya line from the entire line profile. In Paper I we used a single lated Lya profiles for phase 0.26 are compared in Figure 2 

Gaussian PSF with an FWHM = 3.57 km s“ ^ corresponding through the use of the interstellar parameters listed in Table 2. 
to 1.20 diodes at Lya. We find that the resulting changes in the The simulated and observed profiles are in excellent agree- 

line width parameter b for all the lines in our reanalysis of the ment, but the LISM parameters differ slightly from those 

phase 0.26 data set lie within the errors of our initial analysis obtained in the original analysis of the phase 0.26 data. We 
(Paper I) and our reanalysis. The parameter b has its usual then analyzed the lower resolution phase 0.80 Lya profiles 

definition b^ = (2/cT/fn) + where ^ is the most probable obtained with the G160M grating. The resulting parameters 

speed of the turbulent motions and m is the atomic mass. for the H i and D i column densities and interstellar param- 

Figure 8 in Linsky et al. (1995) shows the echelle spectra of eters are in excellent agreement with those obtained from the 

the Mg II AA2796, 2803 lines at the opposite orbital phases and higher resolution echelle spectra at phase 0.26, except for the 

the intrinsic emission-line profiles derived in our analysis. temperature which may be a result of inadequate resolution of 

Because the interstellar absorption features are narrow and the the D i line. This gives us confidence that the lower resolution 
S/N of the data is very high, we used a polynomial inter- G160M spectrum did not seriously degrade the quality of the 

polation between the emission-line peaks on either side of the final results. We believe that the small differences in the inter- 

interstellar absorption to produce a smooth and, we believe, stellar parameters between the original analysis of the phase 

very accurate “continuum” against which to measure the 0.26 echelle spectrum and the present reanalysis are a result of 

interstellar absorption. Figure 3 compares the residual inten- both the different PSF used and the small differences in the 

sides at phase 0.80 for the Mg ri lines, obtained by dividing the assumed intrinsic stellar emission lines, 

observed line profiles by the sum of the intrinsic profiles of the Our analysis of the phase 0.80 data and reanalysis of the 
two stars; the computed residual intensities using the inter- phase 0.26 data are in excellent agreement with the same tern- 
stellar parameters are listed in the third column of Table 2. (All perature for the LISM {T = 7000 ± 500 K compared to the 

random errors determined in this paper are ± 2 <r.) The Mg ii previous result of 7000 ± 200 K) and a similar value for the 

column densities, velocities, and broadening parameters turbulent velocity ((J = 1.6 ± 0.4 km s"^ compared to the pre- 
derived from the data at phase 0,80 and rederived from the vious result of ^ = 1.66 ± 0.03 km s“^). We estimate that 

data at phase 0.26 are very close to those determined in Paper I (D/H)lism = 1.60 ± 0.09 x 10“^ for the Capella line of sight, 

for the observations at phase 0.26. where the uncertainty is taken to be the difference between the 



Velocity (km s *) Velocity ^km s *) 

Fig. 3. — Capella at <^ = 0.80 for Mg n k and Mg ii h. Profiles of the interstellar Mg ii k (left) and h {right) lines (dots with ± I a error bars) at = 0.80 after 
dividing by the sum of the assumed intrinsic line profiles of the two stars. Solid lines are simulated interstellar line profiles for the parameters listed in Table 2. 
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TABLE 2 

Interstellar Parameters for Capella and Procyon Lines of Sight 

Capella Procyon* 


Parameter phase 0.26 phase 0.80 broadened O 

Mgii fe y, (km s“M 23.04 + 0.10 20.93 + 0.10 21.53 + 0.10 21.53 + 0.10 

Mg II k 6 (km 2.64 ±0.13 uncertain 3.13 + 0.16 3.13 + 0.16 

Mg II k (10^^ cm”^) 6.44 ± 0.32 uncertain 3.54 + 0.16 3.54 + 0.16 

Mgii /i u, (km s"‘) 22.72 + 0.10 21.27 + 0.10 21.18 + 0.10 21.18 + 0.10 

Mgiik6(kms“') 2.43 ±0.12 2.51 ± 0.12 2.92 ± 0.15 2.92 ±0.15 

Mgii k N^^(10'2 cm~2) 7.03 ±0.35 6.94 ± 0.35 3.23 ± 0.17 3.23 ±0.17 

Fe II (km s" *) 20.65 ±0.10 ... 20.23 ± 0.10 20.23 ±0.10 

Fell 6 (km s‘) 2.36 ±0.12 ... 2.71 ± 0.13 2.71 ±0.13 

Fe II (10*^ cm'^) 3.12 ±0.16 ... 1.86 ±0.10 1.86 ± 0.10 

Hiu, (kms ‘ ) 21.74 ±0.10 22.48 ±0.10 22.00 ± 0.10 22.00 + 0.10 

Hi6(kms“‘) 11.15 ±0.10 10,59 ±0.15 10.83 ±0.15 10.95 ±0.10 

H I Ni,(10^® cm”^) 1.74±0.02 1.73 ± 0.02 1.18 ± 0.03 1.14±0.02 

Dir, (km s"') 22.04 ±0.10 21.48 ±0.10 22.00 ± 0.10 22.00 ±0.10 

Di6(kms"‘) 7.96 ±0.10 7.58 ± 0.10 7.82 ± 0.10 7.91 ±0.10 

DiNi.(10‘^cm-^) 2.85±0.10 2.71 ± 0.10 1.65 ± 0.10 1.82±0.10 

T (K) 7320 ± 120 6567 ± 120 6769 ± 35 6931 ± 110 

<^(kms •) 1.47 ±0.31 1.65 ± 0.20 2.15 ±0.12 2.13 ±0.13 

D/H ( X 10 - ""1 1 .64 ± 0.05 1 .57 ± 0.05 1 .40 ± 0.20 ( 1 .60) 


• The Mg II and Fe ii parameters for Procyon assume a one-component interstellar medium. 

two different determinations. These cited errors are random, detailed simulations and contour plots that demonstrate the 

since they represent the range in the parameters consistent tradeoff between thermal and nonthermal broadening that is 

with the noise of the data. The ±2 o uncertainties in /Vh i and consistent with the uncertainties in the measured b values. 

Z>Hi listed in columns (2) and (3) of Table 2 are derived from There are also systematic errors associated with the uncer- 

contour plots (e.g., Fig. 4) corresponding to different values of tain PSF and intrinsic stellar line profiles. We estimate these to 

between the simulated and observed profiles. The increase in be ±400 K for the temperature and ±0.2 km s"^ for the 

the errors for T and i compared to Paper I results from more turbulent velocity values based on the differences between the 



1.710 1,7^^0 1.730 1,740 1.750 

N(H 1) (10'® cm 

Fig. 4. — Contour plot of the H i column density ^ vs. hydrogen-broadening parameter , to fit the Capella data at 0 = 0.80 
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TABLE 3 


Local Inteeistellar Medium Properties Toward Procyon and Capella 


1 


«HI 

D/H 


(10*’ cm'^) 

(km s 

(cm 

(io‘h 

References 



Procyon (/ = 214", 

,6 - 13", d 

= 3.5 pc) 

10-14 

10-15 

0.09-0.13 

o 

1. 

Anderson et al. 1978 

10-22 

3-15 

0.09-0.20 

>0.8 

Murthy et al. 1987 

11.2-11.8 

10.68-10.88 

0.104-0.109 

(1.6) 

This paper 



Capella (/ = 163", 

h = 5\d=^ 

12.5 pc) 

9-15 

10 

0.02-0.04 

22-9.6 

Dupree, Baliunas, & Shipman 1977 

15-19 

8-10 

0.04-0.05 

1. 8-4.0 

McClintock et al. 1978 

9-17 

3-12 

0.022-0.04 

^2.0 

Anderson 1979 

6-18 

13-17 

0.014-0.046 

2.4-6.0 

Murthy et al. 1990 

17-21 

10.9-11.4 

0.044-0.055 

1.47-1.72 

Linsky et al. 1993 

17.1-17.6 

10.44-11.25 

0.0443-0.0456 

1.41-1.74 

This paper 


derived interstellar parameters for the two different PSFs. The 
uncertainty in D/H is smaller than that given in Paper I, since 
we now have data at two different phases with which to infer 
the intrinsic line profiles. We estimate the systematic error in 
(D/H)lism to be ( + 0.05, —0.10) x 10“^. Table 3 compares 
these interstellar gas properties with previous determinations 
for the Capella and Procyon lines of sight, 

3. THE PROCYON LINE OF SIGHT 
3. 1 . Observations of Procyon 

Our second target was Procyon, an F5 IV-V star located 3.5 
pc along a line of sight about 52° from Capella. We observed 


the target on 1992 December 21 in the same way as we 
observed Capella at phase 0.26, except that the Lya line was 
observed through the small science aperture (SSA) with the 
G160M grating on side 2. The spectral resolution at Lya was 
only 20,000 (15 km s”^) compared with 84,000 (3.57 km s"^) 
when we used the Echelle-A grating. The Mg ii and Fe ii lines 
were observed through the SSA with Echelle-B, providing the 
same high spectral resolution as for Capella, These obser- 
vations were obtained and reduced in the same way as the 
phase 0.80 Capella observations, except that we used the 1993 
September version of CALHRS. 

Figure 5 illustrates the echelle spectra of the Procyon Mg ii 
resonance lines. The broad absorption line wings are produced 
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Fig. 5. — Echelle-B spectrum of Procyon including the Mg ii lines. Narrow absorption lines to the right of the centers of each self-reversed emission line are 
interstellar. 
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by the optically thick Mg ii lines in Procyon’s photosphere. 
The bright emission features centered at 2796 and 2803 A are 
produced by collisionally excited Mg ii in the warm chromo- 
sphere. The absorption features in the centers of the emission 
lines are self-reversals produced by Mg ii ions at the top of the 
chromosphere. Here the populations of the excited states are 
smaller than at the base of the chromosphere because of 
decreased collisional excitation at the lower densities and the 
lower line radiation field that occurs where the line becomes 
effectively thin. About -1-0.2 A from the central absorption 
features of both lines are very narrow absorption features pro- 
duced by Mg II ions in the interstellar medium. Figure 6 shows 
the photospheric Fe ii A2599 absorption line upon which a 
narrow interstellar component also is superposed. 

3.2. Analysis of the Procyon Line of Sight 

Figures 7 and 8 illustrate the Mg ii and Fe ii residual inten- 
sities obtained by dividing the observed absorption lines by an 
interpolated stellar emission-line profile. We initially assumed 
a single velocity component for the Procyon line of sight and 
deduced the interstellar Mg ii and F ii column densities, 
velocities, and broadening parameters that minimized the 
residuals of the fit. Figures 7 and 8 compare the observed and 
simulated absorption profiles for the one-component model of 
the LISM that minimizes of the residuals. The parameters 
for this model are given in column (4) of Table 2. 

Careful inspection of Figure 7 shows good agreement 
between the observed and computed residual intensities on the 
long-wavelength side of all the lines but small systematic dis- 
agreements on the short-wavelength side of the lines. The best- 


fit one-component model predicts Mg ii residual fluxes that are 
too low at 20-21 km s~^ and too high at 14-17 km s"^ As 
shown in Figure 9, we are able to improve the fits to the data 
greatly (the reduced x^ decreases from 3.2 to 1,3) with a two- 
component model in which the bulk velocities differ by 2.6 km 
s"^ (mean value for the h and k lines), slightly less than the 
instrument resolution of 3.54 km s"^ The parameters for the 
best-fit two-component model are given in Table 4. About 64% 
of the Mg II column density is in the component at lower radial 
velocity. The Fe ii line data are too noisy to discriminate 
between the one- and two-component models (see Fig. 8). 

Since the interstellar H i absorption is saturated and broad, 
we must develop a detailed model for the intrinsic stellar Lya 
line profile against which to measure the interstellar H i and 
D I absorption. Fortunately Procyon (unlike Capella) is 
similar to the Sun in its spectral type, low rotational velocity, 
and weak chromospheric emission. Thus, the intrinsic chromo- 
spheric line shape should be similar to that of the Sun. Since 
the Lya lines of solar-type stars broaden with increasing lumi- 
nosity (Landsman & Simon 1993), we broadened the solar line 
profile (about line center) by a range of multiplicative factors to 
best match the observed line profile outside of the interstellar 
absorption. A broadening factor of 1.4, as shown in Figure 10, 
produces a good fit to the observed profile, except in the region 
of the D I line. 

We have fitted the H i and D i lines with both one- 
component and two-component models for which the veloci- 
ties and relative column densities of each component were 
assumed to be the same as for Mg ii. The resulting interstellar 
parameters and H i and D i column densities of the one- 
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Fig. 6. — Echelle-B spectrum of Procyon including two Fe ii resonance lines. Narrow absorption line at 2599.6 A is interstellar. 
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Fig. 7. — Echelle-B spectra of Procyon’s normalized interstellar Mg ii k (2796 A) {left) and h (2803 A) {right) lines (with ± 1 (7 error bars): the profile that best fils 
the data using the single velocity component interstellar parameters in Table 2 {solid line), and the best-fit profile before instrumental smoothing {dashed line). 
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Fig. 8. — Echelle-B spectrum of Procyon’s normalized interstellar Fe ii A2599 line (with ± 1 o error bars). Left, best fit to the data using the one velocity 
component interstellar parameters in Table 2 {solid line), and the best fit before instrumental smoothing {dashed line). Right, same as left panel, but for the two velocity 
component parameters in Table 3. 
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Fig. 9.~EchelIe-B spectra of Procyon’s normalized interstellar Mg ii k (2796 A) (left) and h (2803 A) (right) lines (with ± I <r error bars), the best fit to the data 
using the two velocity component interstellar parameters in Table 3 [solid line), and the best fit before instrumental smoothing (dosfied line). 


component models (col. [4] in Table 2) and two-component 
models (Table 4) are nearly identical, except that the b values 
for the two-component models are smaller, especially for the 
Mg II and Fe ii lines, and the derived value of ^ is smaller. 
These results are expected because the instrumental resolution 
and the widths of the H i and D i lines are much larger than 
the velocity separation of the two interstellar components, but 
for Mg II and Fe ii the instrumental resolution and line widths 
are similar to the interstellar velocity separation. 

TABLE 4 


Two-Component Interstellar Parameters for the Procyon 
Line of Sight 


Parameter 

Component 1 

Component 2 

Mg 11 k (km s“ b 

21.24 ± 0.10 

24.04 ± 0.10 

Mg II k 6 (km s “ b 

2.3 ± 0.10 

2.3 ± 0.10 

Mg II k Nf. (10^^ cm-2) 

2.29 ±0.10 

1.4 ± 0.08 

Mg 11 h (km M 

20.35 ±0.10 

22.7 ± 0.10 

Mg ii h b (km 

2.30 ±0.10 

2.30 ± 0.10 

Mg II h (10‘^ cm"^) 

2.29 ±0.10 

1.2 ±0.08 

Fe II (km s" M 

19.0 ± 0.10 

22.0 ±0.10 

Fe II b (km s' ') 

2.1 ±0.10 

2.1 ±0.10 

Fe II cm'^) 

1.12 ± 0.06 

0.88 + 0.05 

H I V, (km s'*) 

20.5 ±0.10 

24.0 ± 0.10 

H I 6(km s'*) 

10.78 ±0.10 

10.78 ± 0.10 

H I N^(10*® cm'2) 

0.75 + 0.02 

0.40 ± 0.01 

D I L’f (km s" *) 

20.5 ± 0.10 

24.0 ± 0.10 

D I b (km s' *) 

7.59 + 0.10 

7.59 ± 0.10 

D 1 N, (10*^ cm'^) 

1.19 ±0.10 

0.64 ± 0.04 

r (K) 

6900 + 80 

6900 ± 80 

{ (km «■') 

1.21 ±0.27 

1.21 ±0.27 


We present in column (4) of Table 2 the H i and D i param- 
eters for the best one-component model. With the assumed 
broadened solar Lya line profile, our fit to the shapes of the H i 
and D i line profiles leads to a hydrogen broadening parameter 
6 hi = 10.83 ±0.15 km which is the same as the average 
value 6 hi = 10.87 ± 0.28 km s" ^ that we now find for Capella. 
The thermal (T = 6769 ± 35 K) contribution to this broaden- 
ing is similar to the rederived values for the Capella line of 
sight, but the turbulent contribution (<J = 2.15 ± 0.12 km s 
is larger. 

We have explored a range of parameters j and D/H to 
minimize which characterizes the difference between the 
observed and simulated Lya profiles. Figure 1 1 shows model 
fits to the observed D i line assuming the broadened solar 
emission-line profile and (D/H)hsm = 1-^’ l-^» 1.6 x 10 

None of the fits are particularly good, although they suggest 
that (D/H)lism lies in the range 1 .40 ± 0.20 x 1 0 “ ^ 

Since our assumption that the intrinsic Procyon Lya line is a 
broadened solar profile did not lead to an acceptable fit to the 
observed D i line for any value of D/H, we instead assumed 
that the interstellar D/H ratio is the same as for the Capella 
line of sight [(D/H)usm = L60 x 10 '^] and modified the shape 
of the broadened solar line profile to minimize between the 
observed line profile and the intrinsic line profile folded 
through the interstellar medium and the instrument. Figure 12 
illustrates our derived intrinsic Lya profile and compares the 
observed and simulated line profiles. In Figure 13 we compare 
the observed line profile to models for(D/H)L,sM = L4, 1.6, and 
1.8 X 10’ ^ Now the match to the observed profile is much 
improved. 




1214,5 1215.0 1215.5 1216.0 1216 

Wavelength (A) 

Pjg 10 . — Solid line: Observed Lya spectrum of Procyon. Long dashes: Solar profile broadened about line center by a factor of 1.4. Short dashes: Broadened solar 
profile folded through the interstellar medium with the parameters in the fourth column of Table 2 including (D/H)l, 5 „ = 1.40 x 10" 
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Fig. 11. — Expanded view of the deuterium line in Procyon's spectrum. Solid line: Observed spectrum. Dashed lines : Broadened solar profile folded through the 
interstellar medium with the parameters in the fourth column of Table 2 and {top to hoftom) (D/H)l,sm - 1.2, 1.4, and 1.6 x 10“ 
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Fig. \2— Solid line: Observed Lya spectrum of Procyon. Long dashes: Solar profile broadened about line center by a factor of 1.4 and modified by enhancing 
the blue peak. Short dashes: Modified solar profile folded through the two velocity component interstellar medium with the parameters in Table 3 and(D/H), nu = 

1.60 X 10-^ 
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Fig. 13. — Expanded view of the D i line in Procyon 's spectrum. Solid line : Observed spectrum. Dashed line : Modified solar profile folded through the two velocity 

componentinterstellarmedium with the parameters in Table 3 and (from lop to bottom) (D/H )lism = l A 16 , and 1.8 x 10 ~’. 
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We consider now whether the asymmetric Lya profile with a 
brighter short-wavelength peak is confirmed by any other evi- 
dence. Neff et al. (1990) found that the Mgn and Lya profiles 
for the G8 in star 6 Lep have the same shape when compared 
on a common AA, not velocity, scale. Ayres et al. (1993) found a 
similar scaling for the Mg ii and Lya profiles of the two 
Capella stars and offered a qualitative explanation for the 
curious effect based on the work of Gay ley (1992). We compare 
in Figure 14 the modified broadened solar Lya line with the 
observed Mg ii k line plotted on the same AA scale. The 
remarkable agreement between the two profiles, especially in 
the region of the D i line, provides evidence that our assumed 
intrinsic Lya line profile may be correct. We conclude, there- 
fore, that the Procyon data are consistent with the Capella 
line-of-sight value (D/H)l,sm = 1.60 x 10 Nevertheless, the 
uncertainty in the intrinsic Lya profile of Procyon prevents us 
from firmly concluding that (D/H)l,sm is constant in the LISM. 

The interstellar gas parameters that we derived by using the 
modified solar profile are listed in the right-hand column of 
Table 2. They lie within the random errors of the values 
derived using the broadened solar profile (col. [4] in Table 2) 
and, except for the values of b and they are similar to the 
values for the two-component model in Table 4. The ±2 a 
random error uncertainties in N^^l and listed in column (5) 
of Table 2 are derived from contour plots (e.g., Fig. 15) corre- 
sponding to different values of between the observed profile 
and those simulated with different interstellar parameters but 
the same assumed stellar emission-line profile. It is difficult to 


estimate the systematic errors in the interstellar parameters 
that arise from the uncertain intrinsic Lya profile and the 
uncertain PSF, but the difference between the values of T 
derived using the two different Lya profiles suggests a system- 
atic error of ± 300 K. We therefore adopt for the Procyon line 
of sight the parameters T = 6900 + 80 ( ± 300 systematic 
error) K and = 1.21 ± 0.27 km s‘L 

4. DISCUSSION 

4.1. Does D/H Vary in the Local Interstellar Medium? 

Our analysis of the high-quality GHRS spectra of the 
Capella and Procyon lines of sight are consistent with a single 
value of (D/H)l,sm = 1 -60 ± 0.09 (-1-0.05, —0.10 systematic 
error) x 10"^ for the local interstellar medium, in agreement 
with the mean value (1.5 ± 0.2 x 10'^) found by McCullough 
(1992). There are as yet no other published values of D/H 
based on GHRS spectra for other lines of sight, although 
several nearby stars have been observed. We will be extending 
our study of D/H to other lines of sight in upcoming GHRS 
observations. Eventually, the Far-Ultraviolet Spectroscopic 
Explorer (FUSE) will extend the program to the higher Lyman 
lines, and thus to larger H i columns and more distant targets. 
If real variations in D/H are found, then we will have a power- 
ful tool for studying the physical processes that drive the 
chemical evolution of the Galaxy over short distances and 
timescales. Bruston et al. (1981) have discussed mechanisms 
that could explain variations in D/H on short distance scales. 
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FtG. 14. — Histogram: Observed Lya profile of Procyon. Long dashes: Modified solar profile. Short dashes: Observed Procyon Mg n k line plotted on the Lya 
wavelength scale and scaled to the same peak flux as the modified solar profile. 
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Fig. 15. — Contour plot of H i column density vs. hydrogen-broadening parameter 6 h j to fit the Procyon data by using the modified solar profile. 


4.2. Properties of the Local Interstellar Gas 

Analyses of optical and ultraviolet absorption lines toward 
nearby stars (including late-type stars, hot white dwarfs, and a 
few B-type stars), studies of interstellar gas flowing into the 
heliosphere, and recent studies of EUV absorption toward 
nearby stars support the currently accepted model for the 
LISM within the region of space where < 10^^ cm“^. In 
this model the Sun is located within, but near the edge of, a 
warm cloud (also called the local cloud or the local ‘‘fluff”) in 
which the gas is mostly neutral with «h i ^ 1 cm ~ ^ and 

T 10^ K (see, e.g., Frisch & York 1983; Bruhweiler & Vidal- 
Madjar 1987; Frisch 1994a, b; Vallerga et al. 1994). The local 
cloud extends for roughly 60 pc (to ^ 2 x 10^^ cm'^) 
toward the Sco-Cen Association but only 2-3 pc in some other 
directions. The minimum value of /Vhi 1*cs in the range 0.5- 
1.0 X 10‘® cm “ ^ (Vennes et al. 1994). The local cloud is embed- 
ded in hot ionized gas (T ^ 10^ K) with much lower density 
(«Hi ^ 10“^'^ cm^^) and separated from it by a thin conduc- 
tive interface. Other warm clouds are seen in some directions. 
Thus, lines of sight toward some stars located outside the local 
cloud show no significant additional H i column beyond that 
from the local cloud itself. 

In their analysis of high-resolution Ca ii K line absorption 
toward six nearby stars, Lallement & Bertin (1992) derived the 
flow vector for the local interstellar gas. They found a primary 
bulk flow velocity of 25.7 km s~ ^ from the direction of Galactic 
coordinates I = 186? I and h = — 16?4. Two additional absorp- 
tion components are present along the 5 pc line of sight toward 
a Aql (Ferlet, Lallement, & Vidal-Madjar 1986), and one addi- 
tional component is seen along the 2.7 pc line of sight toward 
Sirius (Lallement et al. 1994). The Capella spectra show only 


one velocity component at 22.0 ± 0.9 (heliocentric) km s ^ 
which agrees precisely with the 22.0 km s"^ projection of the 
local flow velocity vector toward Capella. Both a Aql and 
Sirius also show an absorption component at the projected 
velocity of the local cloud. For Procyon we find two velocity 
components: component 1 at 20.8+ 1.5 km s“^ and com- 
ponent 2 at 23.4 ±1.5 km s“^ based on the means of the 
measured velocities of the Mg ii h and k lines. The 1 a error of 
+ 1.5 km s“^ is the expected absolute velocity error for GHRS 
echelle spectra obtained through the SSA with platinum lamp 
calibration spectra (see Paper I). The velocity errors in Tables 2 
and 4 are random measurement errors that do not include 
systematic effects. We do not use the Fe ii data which are much 
noisier than the Mg ii data. Component 1 is consistent with the 
19.8 km s~ ^ projection of the local cloud flow vector found by 
Lallement & Bertin (1992). Component 2 is not from the same 
cloud as the second component found by Lallement et al. 
(1994) toward Sirius, which is in the same portion of the sky as 
Procyon, since component 2 is +2.6 km s“‘ relative to the 
local cloud, whereas the second component toward Sirius is 
— 6.2 km s“ ^ relative to the local cloud. 

The kinetic temperature of the warm gas in the local cloud 
has been estimated by several techniques. Recent in situ mea- 
surements of interstellar gas flowing through the heliosphere 
yield temperatures of 7000 + 1000 K for helium and 
8000 + 1000 K for H I atoms (see discussion in Lallement & 
Bertin 1992 and Paper I). Measurements of absorption-line 
widths can lead to overestimates of the gas kinetic temperature 
if unresolved velocity components or turbulence are present. 
This may explain why Gry, York, & Vidal-Madjar (1985) 
reported T = 11,000-12,500 K for the gas along the line of 
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sight toward p CMa, which is about 200 pc distant, from their 
analysis of Copernicus Lyman line absorption spectra. 

We believe that our GHRS spectra provide the most accu- 
rate temperatures to date for the warm neutral gas of the 
LISM, because with a spectral resolution of 3.54 and 3.27 km 
s“^ for the high-opacity Mg ii and Fe ii resonance lines, 
respectively, we can find only one velocity component for the 
Capella line of sight and two components toward Procyon. 
These spectral lines of high-mass ions accurately measure the 
nonthermal (turbulent) broadening. At the same time, the 
width of the D i Lya line and the shape of the H i absorption 
close to the saturated core tightly constrain the thermal 
broadening, especially for the echelle spectrum of Cap>ella. The 
agreement between the temperature and turbulent velocities 
for the Capella line of sight, as inferred separately from the 
echelle and G160M spectra of the Lya lines (see Table 2), indi- 
cates that the lower spectral resolution of the G160M data 
does not greatly degrade the accuracy of these results. We find 
that the gas along the Capella and Procyon lines of sight has 
the same temperature (about 7000 K), but the measured turbu- 
lent velocity for the Procyon line of sight is smaller, most likely 
due to the identification of two bulk velocity components 
along the line of sight. We propose, therefore, that the Procyon 
value of ^ = 1.21 ± 0.27 km s“ ^ is the best choice for the turb- 
ulence in the local cloud. This value is consistent with ^ = 1.4 
(-1-0.6, — 1.4) km s' ^ that Lallement et al. (1994) found for the 
local cloud gas toward Sirius and is consistent with the values 
T = 7200 ± 2000 K and ^ = 2.0 ± 0.3 km s“ ^ that Gry et al. 
(1995) recently obtained for component 1 in the line of sight 
toward e CMa. 

We find that the mean H i density is a factor of 2.4 larger for 
the Procyon line of sight («hi = 0.1065 ± 0.0028 cm~^) com- 
pared to the Capella line of sight («h i = 0.0450 ± 0.0006 
cm"^). Since the directions toward Procyon and Capella are 
not very far apart, we can simply explain the different densities 
by placing Procyon near the edge of the local cloud and 
Capella well outside the local cloud. If this geometry is valid, 
then the derived density toward Procyon provides an accurate 
measurement of the local cloud density, but the cloud extends 
only 5.3 pc toward Capella with the remaining 7.2 pc of the line 
of sight consisting of hot gas with essentially no neutral H i. 
The existence of two clouds along the line of sight to Procyon 
makes our determination of the local cloud density and the 
extension of the local cloud toward Capella uncertain. 

4.3. The Boundary Layer between the Local Cloud 
and the Hot Interstellar Medium 

Why do Procyon and Sirius, two very nearby stars separated 
in the sky by only 25'', show second velocity components that 
differ from each other in radial velocity by 10.4 km s“\ while 
Capella, which is separated from Procyon by 52°, does not 
show a second velocity component? If we make the sensible 
assumption that a second velocity component is a result of 
absorption by gas in a distinct cloud along the line of sight to 
the star but outside of the local cloud, then we can place con- 
straints on the cloud size. If we place the Sirius cloud close to 
the 2.7 pc distance of Sirius, then its dimension toward the 
Procyon line of sight must be less than 1.2 pc, or we would see 
an absorption feature at the appropriate velocity in Procyon’s 
spectrum. Similarly, the Procyon cloud cannot extend more 
than 1.2 pc in the direction of Sirius, and both clouds cannot 
extend to the Capella line of sight which is somewhat more 
distant. This implies that both clouds have dimensions that are 


roughly 1 pc or smaller and that they lie at or just beyond the 
edge of the local cloud. 

We believe that the GHRS spectra of the three stars provide 
us with the first glimpse of the boundary layer between the 
local cloud and the surrounding hot interstellar gas. We specu- 
late that this boundary layer differs greatly from the usual 
picture of a smooth conductive interface. Instead, the bound- 
ary layer appears to be turbulent, consisting of many small 
cloudlets that are being sheared from the local cloud by shocks 
or uneven ram pressure as the local cloud and the hot inter- 
stellar gas collide. These cloudlets presumably have conductive 
interfaces with the hot gas and may be destroyed on short 
timescales due to this heating and shearing to even smaller 
scales. Additional observations of other nearby stars are 
needed to confirm this cloudlet boundary layer hypothesis. 

4.4. A Speculation Concerning the Origin 
of the Nonthermal Motions 

Until now, we have characterized the nonthermal motions in 
the local cloud as Gaussian “turbulence” that adds quadrati- 
cally with the thermal motions to explain the observed line 
broadening. The magnitude of the turbulence is very subsonic 
which is a little surprising, since it implies that the thermal and 
turbulent motions are not equilibrated and there are no shock 
waves in the local cloud. Shocks can be produced easily by 
cloud-cloud collisions. We consider, therefore, an alternative 
explanation for the nonthermal line broadening. 

The very high S/N of our Mg ii line spectra provided the 
opportunity to discover two velocity components along the 
Procyon line of sight that diflfer by 2.6 km s“\ close to 
the instrumental resolution of 3.54 km s~\ Similar quality 
data provided no indication of a second velocity component 
along the Capella line of sight, although additional com- 
ponents at smaller velocity separations are possible. Our cloud 
model with a single bulk velocity vector is, of course, an ideal- 
ization. The local cloud may not be rigid, but instead may have 
a small velocity gradient (or velocity shear) such that over 
distances of a few parsecs the line-of-sight velocity changes 
gradually by a few km s ^ This velocity gradient may be 
concentrated near the cloud boundary in certain directions. 
Support for this speculation can be found in the identification 
by Lallement et al. (1994) of cloud G, located very close to the 
local cloud with a very similar flow direction and a flow speed 
that differs by only 3 km s“ ^ from that of the local cloud. They 
argue that the Sun has already moved most of the way through 
the local cloud and will next enter cloud G. We suggest that the 
local cloud and cloud G are not isolated entities, but rather 
parts of the same cloud that contains small velocity gradients. 
We speculate that the measured nonthermal velocities, which 
we have called “turbulent,” may instead indicate the magni- 
tude of the velocity gradients along the line of sight. 

4.5. Gas-Phase Abundance of O in the Capella Line of Sight 

In addition to the interstellar H i, D i, Fe ii, and Mg ii 
absorption lines discussed in this paper and in Paper I, Linsky 
et al. (1995) also found interstellar O i Al 302.2 absorption in 
Capella’s spectrum. Since this feature was observed with the 
G140M grating through the large science aperture (LSA) (see 
Fig. 4a in Linsky et al. 1995), the spectral resolution is insuffi- 
cient for the type of detailed spectral synthesis that we used to 
determine the column densities and broadening parameters for 
the H I, D I, Fe ii, and Mg ii lines. Instead, we infer the inter- 
stellar oxygen abundance from the measured equivalent width 
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of the O I A 1302.2 absorption line, = 0.052 ± 0.002 A. The 
error in Wqi is very small because the data have very high S/N 
and the stellar O i emission line, which forms the “ continuum ” 
against which the absorption is measured, varies smoothly on 
either side of the narrow interstellar feature. 

Since the spectral resolution does not permit us to measure 
the O I Doppler parameter directly, we infer ho i = 3. 1 3 km s " ^ 
from the temperature and turbulent velocity for the Capella 
line of sight ( T = 7000 K and (^ = 1 .6 km s ' ' ) (see Table 2). We 
use /= 0.04887 (Morton 1991) for the oscillator absorption 
strength of the O i line and the curve-of-growth techniques 
outlined by Spitzer (1978) to calculate an O i column density 
of Nqi = 8.2 X 10^^ cm"^ O i is strongly coupled to H i via 
charge exchange collisions. Therefore, O i and H i will exist in 
the same regions of space, and Nqi/N^] should be a good 
measure of the gas phase abundance of oxygen. The oxygen 
abundance is then log (No//Vh) = —3.32, assuming an Hi 
column density of Nhi = 1 735 ± 0.025 x 10^® cm ^ (see 
Table 2). The solar abundance is log (0 /H)q = —3.07 (Anders 
& Grevesse 1989), so our measured O abundance represents a 
logarithmic depletion of D(0) = log (Nq/N„) — log (0 /H)q = 
— 0.25. The O depletion in the LISM toward Capella is similar 
to the measurements for other, more distant lipes of sight 
(Jenkins 1987). We also note that the O abundance toward 
Capella is the same as that measured for B star photospheres 
(Gies & Lambert 1992). The issue of choosing the “cosmic*’ 
reference abundances for interstellar studies is an important 
one (see § 4.2 of Sofia, Cardelli, & Savage 1994). Additional 
observations of O i in the LISM may provide important new 
insights regarding this problem. 

4,6. Gas- Phase Abundances of Fc and Mg in the Procyon 
Line of Sight 

We now estimate the abundances of Fe and Mg in the 
Procyon line of sight by making the same assumptions as in 
Paper I. For the column densities of Fe ii. Mg ii, and H i we 
sum the contributions of the two velocity components: 
Npeii = 2.0 X 10^^ cm■^ Nm.ii = 3.6 x 10^^ cm and 
Nhi = 1.15 ± 0.03 X 10^® cm* , respectively. Therefore, log 
(^Fcii/^Hi) = -5.76 and log (Nm.h/Nhi) = “5.50. The Fe ii 
abundance for the Procyon line of sight is the same as for 
Capella, but the Mg ii abundance is 0.06 dex smaller. Compar- 
ing these numbers to the abundances of Fe and Mg, log 
(Fe/H)o = -4,49 and log (Mg/H)© = -4.41 from Anders & 
Grevesse (1989), and assuming that the contributions to Fe ii 
and Mg ii from ionized gas are negligible, we derive the 
logarithmic depletions, D(Fe) = (A^fcii/^hi) “ (Fe/H)© = 

-1.27 and D(Mg) = (Nm,„/N„,) - log (Mg/H)© = -1.09 in 
the warm gas toward Procyon. 

These Mg and Fe depletion values can be compared to 
Copernicus estimates of £)(Mg) = —0.3 and D(Fe) = —1.4 in 
the warm interstellar medium (Jenkins, Savage, & Spitzer 
1986). The value for Fe is similar to that measured for the 
Capella and Procyon sight lines, while the value for Mg differs 
by 0.8 dex. Most of the difference may be due to a systematic 
error in the/- values for the far-ultraviolet Mg ii lines near 1240 
A used in the Copernicus study. Sofia et al. (1994) found from a 
careful analysis of GHRS echelle mode data for interstellar gas 
toward ^ Per that the far-UV lines of Mg ii yield Mg ii column 
densities 0.67 dex larger than those obtained from an analysis 
of the strong but damped Mg ii lines near 2800 A. Since the 
/-values of the Mg ii A2800 lines are well determined, they 
recommended a 0.67 dex revision to the /-values for the Mg ii 


far-UV doublet. Although the Mg ii results may be uncertain, 
the measurements for Fe ii clearly reveal that substantial gas- 
phase depletion occurs in the cloud toward Capella and 
Procyon. Evidently the Fe in the cloud is mostly found in 
interstellar dust. 

5. COSMOLtXJICAL SIGNIFICANCE OF THESE RESULTS 

As described in § 1, an important objective of our program is 
to measure D/H in the LISM and to infer the primordial 
number ratio (D/H)p, which tells us the number density of 
baryons during the period 100—1000 s after the big bang when 
the universe had cooled enough for the light nuclei to form. 
Since the Hubble expansion relates the baryon densities then 
and now, (D/H)^ also determines the mean baryon density in 
the universe today and the ratio Qg of the local baryon density 
to the critical density (the required density to eventually halt 
the expansion). Thus, (D/H)p is an important parameter for 
experimental cosmology. 

Although our data do not allow us to determine (D/H)^ 
directly, we can infer its value from our measurement of 
(D/H)lism and chemical evolution models for the Galaxy. Steig- 
man & Tosi (1992), for example, have calculated the survival 
fraction of deuterium as the primordial D is converted to 
heavier elements in the cores of stars and the deuterium- 
depleted gas is dispersed into the interstellar medium from 
which later generations of stars are formed. Their calculations 
indicate that (D/H)^ = 1.5-3.0 x (D/H)lism» so that (D/H)^ = 
2.3-5.1 X 10“^. On the other hand, Vangioni-Flam, Olive, & 
Prantzos (1994) have argued that an astration factor (the 
inverse of the survival fraction) of about 5 is more reasonable, 
in which case (D/H)^ % 8 x 10"^. Comparison of the value for 
(D/H)p, assuming the Steigman-Tosi astration factor of 1. 5-3.0 
with recent big bang nucleosynthesis calculations (Walker et 
al. 1991), indicates that j/iq = 3.5-5.8, where 10^® times 
the ratio of nucleons to photons by number. The range in f/jo 
leads to the very important result that 0.05 < Clghlo ^ 0.09, 
where /150 is the Hubble constant in units of 50 km s“ ^ Mpc" ^ 
The Vangioni-Flam et al. astration factor instead leads to 

^ 2.8 and Qghlo % 0.042. Thus, no matter what value one 
assumes for the Hubble constant, Clg < 1, and a universe with a 
cosmological constant of zero consisting only of baryons must 
be very open. Tremaine (1992) and others, however, argued 
that »90% of the universe consists of dark nonbaryonic 
matter on the basis of the dynamical properties of galaxies and 
clusters. Thus, whether Q > 1.0 is true remains an open ques- 
tion. 

Another important result in experimental cosmology is that 
our lower limit on (D/H)l,sm ^ F41 x 10“^ provides a hard 
lower limit to (D/H)^ and thus a hard upper limit on 7.8 
and Qghlo < 0.125. These limits are independent of Galactic 
chemical evolution models and only assume that D is 
destroyed with time. 

6. CONCLUSIONS 

In this second in our series of studies of the D/H ratio and 
other properties of warm interstellar gas along the lines of sight 
toward nearby stars, we analyze GHRS spectra toward 
Procyon (3.5 pc) and Capella (12.5 pc). Observations of Capella 
near opposite orbital quadratures (orbital phases 0.26 and 
0.80) permit us to derive both the intrinsic stellar Lya line 
profiles and the interstellar absorption lines more precisely 
than before. We conclude that (D/H)hsm = 1-bO + 0.09 ( -1-0.05, 
-0.10 systematic error) x 10“^. Our analysis of the Lya, 
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Mg II, and Fe ii lines indicates that the temperature of the 
interstellar gas is T = 7000 ± 500 (±400 systematic error) K 
and the microturbulence is ^ = 1.6 ± 0.4 (±0.2 systematic 
error) km s“ *. 

Procyon’s Mg ii line profiles demonstrate that there are two 
velocity components in that line of sight separated by 2.6 km 
s ^ slightly less than the instrumental resolution. Since we 
could not fit the observed D i line shape by assuming that the 
intrinsic Lya line profile is simply a broadened solar profile for 
any value of D/H, we assumed that the D/H ratio is the same 
as for the Capella line of sight and modified the broadened 
solar profile to achieve excellent agreement between the simu- 
lated and observed line profiles. The resulting asymmetric 
intrinsic stellar line profile is similar to the scaled Mg ii profile 
shape, which supports our choice of this Lya profile. The 
Procyon data thus provide a consistency check on the Capella 
value for (D/H)lism = 1*6 x 10“^, but the uncertainty in the 
intrinsic Lya emission-line profile does not permit us to con- 
clude that the D/H ratio is constant in the LISM. Observations 
of other lines of sight are needed to determine whether a single 
value of D/H characterizes the LISM. 

We find that interstellar gas properties of the lines of sight 
toward Procyon and Capella are very similar. Since the mean 
H I density is a factor of 2.4 larger for the Procyon line of sight 
than for Capella, we propose that the first 5.3 pc along the 
Capella line of sight lies within the local cloud, while the 
remaining 7.2 pc lies in the hot gas surrounding the local cloud. 
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This model leads us to propose that = 0.1065 ± 0.0028 
cm“^ be adopted for the neutral hydrogen density within the 
local cloud and ^ = 1.21 ± 0.27 km s“^ be adopted for the 
nonthermal motions. We speculate, however, that the mea- 
sured nonthermal motions which we call “turbulence'’ are 
instead indicative of systematic velocity shear within the local 
cloud: it is not a rigid comoving structure. We find that the 
depletions of Fe and Mg for both lines of sight are at the lower 
end of the distribution for stars, as is the depletion of O along 
the Capella line of sight. 

Using the Galactic chemical evolution calculations for D by 
Steigman & Tosi (1992) or by Vangioni-Flam et al. (1994), we 
infer that 0.05 < ^ ^ 0.042, respectively. 

Thus, no matter what value one assumes for the Hubble con- 
stant, we conclude that 1 and that a universe with a 
cosmological constant of zero consisting only of baryons must 
be very open. The assumption of no Galactic chemical evolu- 
tion leads to a model-independent upper bound to D^/iso < 
0.125. 
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